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Abstract. We fit the single-hadron transverse-momentum spectra measured in 
Pb-l-Pb collisions at s/snn = 2.76 TeV with the blast-wave model that includes 
production via resonance decays. Common fit to pions, kaons, (anti)protons, and 
lambdas yields centrality dependence of the freeze-out temperature and transverse 
expansion velocity. In the most central collisions we see T = 98 MeV and (vt) = 0.654. 
The K* resonance fits into this picture but the (p meson might freeze-out a little earlier. 
Multistrange baryons seem to decouple at higher temperature and weaker transverse 
flow. Within our model we observe hints of chemical potential for the charged pions. 


PACS numbers: 25.75.-q,25.75.Dw,25.75.Ld 

1. Introduction 

Hot nuclear matter created in ultrarelativistic heavy-ion collisions exhibits its properties 
through the collective expansion of the hrebalh In addition to the longitudinal flow 
which is largely set by the initial transparency of the nuclei, pressure gradients drive 
strong expansion also perpendicularly to the collision axis. Due to expansion and cooling 
the hreball eventually arrives to the state where it disintegrates. This is the moment 
when most abundant hadrons obtain their observed momenta and momentum spectra 
are frozen-out. In fact, it does not happen instantaneously for all hadrons. Freeze-out 
is a continuous process. 

It is useful to describe the freeze-out state of the hreball by hydrodynamically 
inspired parametrisations. They usually assume that freeze-out happens instantaneously 
along some hypersurface. The models include locally thermalized momentum 
distribution and some pattern of longitudinal and transverse expansion how velocity. 
A few parametrisations are on the market [D El El H El i which diher in details of the 
assumptions for the density and velocity prohles or the exact shape of the freeze-out 
hypersurface. In our analysis we shall make use of the so-called blast-wave model [3]. 

An important part of hnal state hadrons does not originate directly from the hreball 
but comes from decays of unstable resonances which themselves were emitted from 
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the fireball. Hadron production by resonance decays is clearly substantial for pions. 
However, it also contributes largely to protons and kaons, as we will show below. In 
spite of this, it is often omitted in calculations and analyses due to its computational 
complexity mnu- 

In this paper we compare Monte-Carlo-generated transverse momentum spectra 
with the data measured in Pb+Pb collisions at ^/sNN = 2.76 TeV by ALICE 
collaboration ra- Our simulation follows the prescription of the blast-wave model. 
The Monte Carlo treatment has been chosen because in this way we are able to include 
also the resonance production of hadrons into the simulation. A similar analysis of 
RHIC spectra can be found in HU. 

This allows us to deduce the values of kinetic freeze-out temperature and transverse 
expansion velocity much more reliably than in hts with only direct thermal production 
included, e.g. [XO] . 

In the next Section we briefly summarize the relevant features of our Monte Carlo 
model—the package DRAGON In Section]^ we study how resonance production of 

hadrons affects pt spectra. As a cross-check of our approach we have also performed an 
analysis identical to and then show how results change if resonances are included. 
Our main results are presented in Section We have htted pt spectra of pions, kaons, 
(anti)protons, and A’s [13] and also looked at S’s and D’s [T3|. From the resulting hreball 
freeze-out parameters we calculated spectra of K* and (j) and compared them with the 
measured data |T5|. We show our values of the parameters for six different centrality 
classes. Multistrange baryons seem to drop out of the systematics that describes other 
hadrons. Conclusions are drawn in Section IH 


2. The model 


Our analysis is performed in framework of the blast-wave model, which is characterised 
by its emission function. Formally, it is the Wigner phase-space density of the source of 
hadrons of type i 


S{x,p) d^x = g., 


rrit cosh{ri — y) 


(2i) 


exp 


p^u^^ - Pi 

T 


+ Si 




X r dr dip — tq) t dr dg . 


( 1 ) 


Because of the dominant expansion in the longitudinal direction one uses longitudinal 
proper time r = \/t‘^ — z‘^ and space-time rapidity g = ^ln((f + z)/{t — z)). Polar 
coordinates r, (p are used in the transverse plane. We use proper quantum statistical 
distributions with Sj = 1 (—1) for fermions (bosons) and Pi is the spin degeneracy. Every 
isospin state is treated separately. This prescription assumes sharp freeze-out along the 
hypersurface r = tq and uniform density distribution within the radius R. This means 
that our freeze-out time does not depend on radial coordinate The expansion of the 

f This is one of the differences between our model and the Cracow single freeze-out model [16] which 
was used in fits to pt spectra measured by ALICE collaboration recently [I21II81II9]. 
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fireball is represented by the velocity field 

= (cosh rjt cosh rj, sinh rjt cos (p, sinh rjt sin (p, cosh rjt sinh rf) (2) 

where the transverse velocity is snch that 

/ r \ ” 

vt = tajih.7]t = r]f \^—j . (3) 

In this relation rjf parametrises transverse flow gradient and n the profile of the 
transverse velocity. The mean transverse velocity is then 

The transverse size R and the freeze-ont proper time tq inflnence total normalizations 
of transverse momentnm spectra. However, in this stndy we ignore those and hence we 
have no sensitivity to these geometric parameters. 

From the emission fnnction, spectrnm of directly prodnced hadrons is obtained as 

(5) 

where the integration rnns over the whole freeze-ont hypersnrface. If one replaces the 
qnantum-statistical distribntion in eq. ([^ by the classical Boltzmann distribntion and 
performs some of the integrations in eq. ([^, one arrives at 

^ (^psinh? 7 t(r)j Ki cosh? 7 t(r) j dr. (6) 

This formnla is rather easy to evalnate and thus it is often used in spectra fitting. 

Resonances are emitted as described by the emission function in eq. ([^ and then 
decay exponentially in time according to their width. All matrix elements for the decays 
are assumed to be constant and thus the decay is determined by the phase-space only. 
In principle, one can derive the emission function that includes hadrons from resonance 
decays [3] but the expression—unlike eq. 0 — is not suitable for massively multiple use 
in the fitting procedure. 

Therefore, we generate stable hadrons and also resonances according to distribution 
S{x,p) and the assumption of chemical equilibrium. Resonances are then let to decay so 
that in the end we look only at stable hadrons. In our simulation we will systematically 
vary the temperature T, transverse flow gradient rjf, and the power n. We construct 
histograms in pt for some number of Monte Carlo events for each set of {T,rif,n) 
parameters which are then compared with measured data and a value of is obtained. 
This procedure is rather CPU and storage-demanding because the generation must be 
done separately for each set of parameters that we want to check. 

The share of hadrons coming from resonance decays in our model is set at the 
chemical freeze-out specified by Tch = 152 MeV and ps = 0 m- Note that shifting 
these values slightly is not expected to cause a big change in the shape of the transverse 
momentum spectra. 

This value of chemical freeze-out temperature was found to reproduce the ratios 
of hadron multiplicities [20] and we keep it in order to satisfy that observation. Since 
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the thermal freeze-out temperature is lower, we actually assume here that the time 
between chemical and thermal freeze-out is so short that the majority of resonances does 
not decay so that their products would rescatter. The opposite extreme scenario (not 
studied here) is that strong interactions regenerate some number of resonances. This 
must happen so that the state of partial chemical equilibrium is established, because 
the hnal observed numbers of stable hadrons is hxed to the values that correspond to 
the chemical freeze-out state. Thus in partial chemical equilibrium the loss of pions due 
to fewer resonances that decay would have to be exactly compensated by the increase 
of the number of direct pions. The real scenario may likely be somewhere in between 
these two extreme models. Another possibility is the scenario with single chemical and 
thermal freeze-out [miiH]. As mentioned above, due to fast transverse expansion at the 
LHC, here we assume only short time between chemical and thermal freeze-out which 
keeps basically all resonances according to the chemical freeze-out temperature. At the 
chemical freeze-out, all hadrons are produced from the same common volume. Then, in 
the short time of further expansion the volume occupied by each species may be slightly 
different due to the difference in their thermal velocities. Same kind of model has been 
applied also in an analysis of nuclear collisions at RHIC energies HU. 

3. Impact of resonance decays on transverse momentum spectra 

Data from Pb-|-Pb collisions at ^/snn = 2.76 TeV were measured by the ALICE 
collaboration. Spectra of the most abundant species—pions, kaons, protons and 
antiprotons—were published in uni along with a simple thermal £t to the spectra which 
included only directly produced particles and no resonance decay effect. An effort has 
been made to identify hducial intervals in pt where the influence of resonance decays 
can be omitted. They were 0.5 GeV/c < pt < 1 GeV/c for pions, 0.2 GeV/c < pt < 
1.5 GeV/c for kaons and 0.3 GeV/c < pt < 3 GeV/c for protons. As we first wanted 
to set the benchmark and compare our method to other analyses, we performed the 
same blast-wave model £t with the same hducial intervals as [10] using our Monte Garlo 
approach instead of eq. 0. 

Event samples were generated for different sets of parameters {T,rif,n). The 
parameter space was sampled so that the step in temperature was 4 MeV, the step 
in rjf was 0.01, and the step in n was 0.02 or 0.04. Each of these steps induces roughly 
the same amount of change in the resulting spectra. The step size was compromised 
in view of the GPU time required and has impact on the precision with which we can 
locate the parameter set that fits the data best. Explored parameter sets were chosen 
so that we always explore large enough region in parameter space around the best fit. 

We checked that the generated spectra are accurately described by eq. (§ with 
appropriate parameters. Then we looked for the parameter set (T, rif,n) for which Monte 
Garlo data show the best match with measured data using the minimum method. 
Results are summarised in Table They agree with those by ALIGE collaboration [TO] 
reasonably well. 
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Table 1. Freeze-out parameters obtained from fits to pion, kaon and proton p* spectra 
within fiducial intervals as defined in m- Results without and with resonance decays 
included are compared to the results by ALICE collaboration. 


centrality 


ALICE [TU] 


no resonances 

with 

resonances 

(MeV) 

n 

T 

(MeV) 

(n) 

n 

T 

(MeV) 

(n) 

n 

0-5% 

95 

0.651 

0.71 

98 

0.645 

0.73 

82 

0.662 

0.69 

5-10% 

97 

0.646 

0.72 

98 

0.645 

0.73 

94 

0.654 

0.69 

10-20% 

99 

0.639 

0.74 

102 

0.637 

0.73 

90 

0.649 

0.71 

20-30% 

101 

0.625 

0.78 

102 

0.624 

0.79 

98 

0.633 

0.75 

30-40% 

106 

0.604 

0.84 

no 

0.605 

0.81 

102 

0.616 

0.79 

40-50% 

112 

0.574 

0.94 

no 

0.572 

0.97 

118 

0.581 

0.89 

50-60% 

118 

0.535 

1.10 

122 

0.527 

1.15 

126 

0.541 

1.03 

60-70% 

129 

0.489 

1.29 

126 

0.484 

1.39 

146 

0.489 

1.23 

70-80% 

139 

0.438 

1.58 

142 

0.439 

1.51 

170 

0.423 

1.55 


Then, in order to check whether resonances can indeed be neglected in the hducial 
range, we added resonance decays to our Monte Carlo data. Their contribution is 
characterised by the chemical freeze-out temperature as described in the previous 
section. We still kept the limitation to the hducial pt intervals. Results are summarised 
in Table [Tj The extracted kinetic freeze-out temperatures changed and the direction in 
which they have gone depends on centrality. For central collisions the temperature was 
lowered while for peripheral ones it went up. It thus seems that leaving out resonances 
is not a good approximation even if one tries to identify a reasonable hducial range for 
the htting. 

To support this statement we have analysed our Monte Carlo events and looked 
at the number of hadrons emitted directly and the number coming from the decays of 
resonances. 

In Fig. we plot the ratio of hadrons produced from the resonance decays to 
those produced directly as functions of pt- They are calculated for model parameters 
which best ht the spectra in most central collisions (see Table [^. The portion of 
resonance production is very important not only for pions but also for kaons and protons. 
Moreover, the ratio shows nontrivial pt dependence. One might omit resonances in the 
htting procedure in a region where these ratios are constant but there is no such interval 
for any of the investigated species. For pions the pt-dependence is even non-monotonic. 
Hence, leaving out resonances is an invalid simplihcation. 

We analysed the resonance contribution to hadron spectra in much detail. In Fig. 
we demonstrate that the relative contribution of resonance decays to pion production 
also depends on kinetic freeze-out temperature and transverse how. Simulated data 
correspond to best hts in the indicated centrality classes (next Section). The rise of 
decays/direct ratio with increasing pt is due to decays of resonances with masses so high 
that pions get large kinetic energy from their decays. In central collisions the freeze- 
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Figure 1. Ratios of hadrons coming from resonance decays to those produced directly 
for three different species. Monte Carlo data were simulated by DRAGON with the 
parameters T = 98 MeV, rjf = 0.88, and n = 0.69. 



Pt [GeV/c] 


Figure 2. Ratios of pions coming from resonance decays to those produced directly. 
Monte Carlo data were simulated by DRAGON with the parameters given by the best 
fits in the indicated centrality class. 


out temperature is low. Thus at higher pt the pion spectrum is dominated by heavy 
resonance decay products. It also helps that these resonances are boosted due to large 
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transverse velocity of the part of fireball which emits them. In peripheral collisions the 
freeze-out temperature increases and the transverse flow is weaker. The former gives 
high pt to direct pions and the latter gives less boost to heavy resonances. As a result, 
the contribution of resonance production decreases with growing pt. 

4. Fits to hadron transverse momentum spectra 

In our actual analysis we fitted spectra of pions, charged kaons, protons and antiprotons 
[In] together with spectra of i7°’s and A’s [I3]. Our simulations were further compared 
with spectra of multistrange baryons S and O published in [T3]. Finally, we compared 
with data [15] our predictions for K* and 0 which were based on the results of fits to 
vr, K, p, p, and A. 

The and hyperon data are divided into fewer centrality classes than data from 
7r±, K^, p, and p. In order to make consistent fits which include all species, we used 
the centrality classes of hyperon spectra and summed up corresponding centrality bins 
in pion, kaon, and proton spectra. Error bars of experimental data points were taken 
as statistical and systematic added in quadrature. 

The parameter space {T,pf,n) was sampled as described in the previous section. 
After having all ensembles of events simulated we compared them with data. Histograms 
of simulated data were normalized individually so that their normalization matches that 
of experimental spectrum. Among all simulated ensembles we found the one which best 
agrees with the data (the minimum x^)- We chose this procedure in order to find the 
best fits to the shape of the spectrum which is sensitive to the dynamical state of the 
fireball at the freeze-out. Their normalization, on the other hand, requires the volume 
which is best reconstructed in common analysis with femtoscopy data. 

Fit results depend on the cuts imposed on the spectra. For each species we specified 
a high-pt cut. This is justihed, because from some value of pt upwards hadrons are 
expected to be more influenced by hard production processes and not described by the 
thermal model we use. To find the cuts we determined ratios of measured-to-simulated 
bin content Ri = for each species i as a function of pt- Then, those bins 

from the high-p* end of the spectra for which the value of Ri was outside the range (0.9, 
1.1) were excluded, effectively introducing high-p^ cuts. (The chosen range was found 
to lead to best fits.) In the next step, Monte Carlo data were again compared with 
the data without excluded bins. This time, new high-pt cuts were determined and the 
procedure was repeated again. The iterations were stopped when they converged to one 
set of parameters and certain set of cuts. 

For pions, we imposed also a low pt cut at 400 MeV/c, thus not including the first 9 
bins into the comparison. We found that the quality of the fit significantly deteriorated 
without this cut. It seems that pion spectra show an enhancement here which might be 
due to non-equilibrium pion chemical potential not included in our simulation. 

Here we did not try to come up with a fit function using Tsallis distribution that 
would cover the whole pt interval EIIE21E3I. Its interpretation in terms of thermal 



T [MeV] 

Figure 3. Positions of best-fit values of T and (vt) for transverse momentum spectra 
of TT^, K^, p, p, , and A are shown by red -l-’s. Results go from central (upper 
left) to non-central (lower right) collisions, with centrality classes indicated in Table 
Around each best fit value we estimate the 99% confidence-level region. By blue x’s 
we show results of fits to spectra of multistrange species, centralities are indicated in 
Table For multistrange species we estimate the 68% confidence-level regions. 


model would be unclear. 

4.I. Spectra of p,p,n^,K^,K^, A 

We first fitted pt spectra for the most abundant species: pious, kaons, (anti)protons, 
7A°’s, and A’s. In Fig. we show positions of the best hts for each centrality. Note 
that in the actual htting procedure we determined three parameters: T, r/y, and n. In 
order to compare results in a relevant way we have plotted here the positions of the best 
hts in just two parameters, T and (ut), with the latter calculated from eq. (|^ for the 
values of n, r// that give the minimum of x^{T,Pf,n). As can also be seen in Table 
with the change of centrality, (vt) sometimes changes due to a change in n and not in 
Pf. In Table the htted values are displayed together with their and the number 
of degrees of freedom. In Fig. they are displayed together with our estimates for 
the uncertainty region for the best ht values. We stress that these are only estimates 
because the evaluation of in discrete points does not allow us to determine them more 
precisely. We observe that the position of minimum determines very precisely the value 
of (vt). For the temperature, the uncertainty is larger. It may range up to 10-15 MeV 
from the minimum. 

We illustrate the quality of the hts in Figs. |4]-[^ For brevity, for the most abundant 
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Table 2. The best parameter values resulting from fits to transverse momentum 
spectra of pions, kaons, (anti)protons, _?f°’s and A’s. 


centrality 

T 

(MeV) 

If 

n 

M 

X^/Ndof 

Ndof 

0-5% 

98 

0.88 

0.69 

0.654 

0.214 

194 

5-10% 

98 

0.88 

0.71 

0.649 

0.266 

197 

10-20% 

106 

0.87 

0.71 

0.642 

0.272 

210 

20-40% 

114 

0.86 

0.81 

0.612 

0.294 

202 

40-60% 

138 

0.82 

0.99 

0.548 

0.347 

195 

60-80% 

166 

0.77 

1.43 

0.449 

0.449 

168 


spectra we only show the comparison with positively charged species and note that 
the comparison looks the same for negative species. For pions (Fig. we clearly see 
the enhancement at low pt- Recall that our htting procedure does not include any 
chemical potential but we account for resonance production. For pions at pt below 
400 MeV/c, resonance contribution grows as pt is lowered. This can be observed in 
Fig.ll Still, although qualitatively this behaviour might agree with what is seen in 
data, quantitatively it is not sufficient to explain them. Resonances tend to populate 
low-pt region but this is not enough to £t the data. We speculate that the solution might 
be in introducing non-equilibrium chemical potential for pions. This would naturally 
occur if the hadron gas chemically freezes out at a temperature around 150-160 MeV 
and then cools down while keeping the effective ratios of individual species constant. 
We estimated that the pion chemical potential at kinetic freeze-out temperature might 
reach values around 100 MeV. This is not enough for Bose-Einstein condensation but 
modihes the spectrum considerably. 

In the higher-pt region the pion spectrum is well reproduced up to about 2 GeV. 
This seems reasonable, as for higher pt we may see signs of hard production. 

Charged kaons (Fig. are well reproduced in a similar pt interval without the 
need to leave out any bins at low p*. We also observe good fits to (anti)proton spectrum 
(Fig.@ stretching out even to about 4 GeV. In general, the agreement becomes slightly 
worse when going away from central collisions. For the most peripheral class charged 
kaons depart from the theoretical curve above 2 GeV/c and protons above 3.5 GeV/c. 

In the £t we also included strange V°’s which were measured up to pt of 4.5 GeV/c 
na. The kaons show (Fig. similar behaviour as their charged isospin partners. When 
comparing the two plots for (Fig. and (Fig. one should notice the different 
Pt intervals in which these spectra are measured. Below the pt cut of 2 GeV also the 
A’s are htted well (Fig. |^. They slightly depart from this agreement above 2 GeV, 
but in a different way than other species: A’s are steeper at high pt than the thermal 
hts. Note that the values of freeze-out parameters were determined in common hts to 
all 8 species, including A’s. We might be seeing here the beginning of the departure of 
strange baryon spectra from the scenario of common freeze-out. Note, however, that 
the departure occurs for pt above 2 GeV/c where thermal hydrodynamically inspired 
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0.0 0.5 1.0 1.5 2.0 2.5 3.0 

Pt [GeV/c] 


Figure 4. Transverse momentum spectra of 7r+ for different centralities. To display 
all spectra in one figure we divide data from non-central collisions by factors 2, 4, 8, 
16, and 32. Lower panels show the experiment-to-Monte-Carlo ratios Ri = 
which demonstrate the agreement of measured and fitted spectra. Different panels 
correspond to different centralities. Horizontal lines indicate the interval 0.9 < Ri < 
1.1 to which we limit our fitting procedure. 


model should not be given too much stress. 
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0.0 0.5 1.0 1.5 2.0 2.5 3.0 


Pt [GeV/c] 


Figure 5. Same as Fig. but for positive kaons. 


4-2. Spectra of. 


and H 


The departure of data from the model fitted in previous subsection is clearly seen when 
we focus on pt spectra of multiply strange baryons: hi, and fl. We show only the 

spectra of baryons in Fig. (F“) and Fig. 10 (hi). Results for antibaryons are similar. 


Note that these species were not included in the previous htting procedure. In order to 
check whether they £t into the same description of the fireball we have simulated their 
spectra in DRAGON with the sets of parameters that come from the fits to vr, K, p, 
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0.0 0.5 1.0 1.5 2.0 2.5 3.0 3.5 4.0 4.5 

Pt [GeV/c] 


Figure 6. Same as Fig. but for protons. 


p, A (Table 1^. These predictions are plotted in Figs, [^and 10 by dashed cnrves. We 
notice that in general the data are steeper than the prediction. The data would thus 
indicate lower transverse expansion velocities and earlier freeze-out of multiply strange 
baryons m- Note, however, that for the two most peripheral classes the agreement of 
data and prediction becomes better. 

In order to quantify this statement we have performed common hts to the four 
multistrange baryon species: hi, and and extracted their values of T, r^/, n 
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0.0 0.5 1.0 1.5 2.0 2.5 3.0 3.5 4.0 4.5 

Pt [GeV/c] 


Figure 7. Same as Fig. but for K°. 


and corresponding (vt). They are listed in Table and compared with the other £t 
resnlts in Fig. This clearly shows the general effect that mnltiply strange baryons 
seem to deconple at higher temperatnre and weaker transverse expansion. The fitted 
spectra are compared to data in Figs. and by solid lines. Dne to higher statistics, fit 
resnlts seem to be driven by the S’s rather than by hl’s. Although it is hard to support 
it statistically, it seems that the hi spectrum is even steeper than the fit results. This 
would indicate even lower expansion velocity at the moment of hi decoupling. 
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0.0 0.5 1.0 1.5 2.0 2.5 3.0 3.5 4.0 4.5 

Pt [GeV/c] 


Figure 8. Same as Fig. but for A. 


Looking at the most peripheral centrality class we notice that the measured data 
seem to be well reproduced by both sets of fit results. Note that the freeze-out 
temperature in peripheral collisions is higher than in the central ones. The spectra 
in those collisions are thus less sensitive to the mass of particles. What we see here is 
an illustration of the ambiguity in determination of T and (vt). Lower temperature can 
be compensated by stronger flow and vice-versa. When spectra depend on the mass, 
this ambiguity is usually resolved by different shapes of pt spectra of different species. 
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0.5 1.0 1.5 2.0 2.5 3.0 3.5 4.0 4.5 


Pt [GeV/c] 


Figure 9. Transverse momentum spectra of for centralities indicated in Table 
To display all spectra in one figure we divide data from non-central collisions by factors 
2, 4, 8, and 16. Dashed curves show predictions based on parameters from Tablej^ For 
the most central bin we show two predictions calculated with parameters for centrality 
classes 0-5% and 5-10%; they differ very little. Lower panels show the experiment-to- 
Monte-Carlo ratios Ri = which demonstrate the agreement of measured 

and fitted spectra. Different panels correspond to different centralities. Horizontal 
lines indicate the interval 0.9 < < 1.1 to which we limit our fitting procedure. 


Figure [^indicates that statistically one can still resolve the ht to multistrange baryons 
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1.0 1.5 2.0 2.5 3.0 3.5 4.0 4.5 


Pt [GeV/c] 


Figure 10. Same as but for ft. 


as the uniquely correct result. 

4-3. Spectra of K* and cf 

ALICE collaboration has also published transverse momentum spectra of resonances 
K* and 0 [15]. Stndying resonances may reveal information abont the dynamics of the 
hreball, particularly its evolution between chemical and thermal freeze-out. In order 
to simulate these resonances we have forbidden them to decay in our Monte Carlo 
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Table 3. The best parameter values resulting from fits to transverse momentum 
spectra of S and 12. 


centrality 

T 

(MeV) 


n 

(vt) 

xVNdof 

Ndof 

0-10% 

126 

0.82 

0.71 

0.605 

0.375 

41 

10-20% 

138 

0.81 

0.85 

0.568 

0.325 

41 

20-40% 

158 

0.78 

0.91 

0.536 

0.417 

37 

40-60% 

174 

0.76 

1.11 

0.489 

0.497 

33 

60-80% 

250 

0.64 

1.35 

0.382 

0.795 

39 


generator since by defanlt their decay is tnrned on. We simnlated their spectra with 
the parameters from Table and compare them with data in Figs. and 1^ Hence, 
they have not been fitted in our treatment. We observe reasonably good agreement of 
theory with data in centrality classes 0-20% and 20-40% and pt below 2 GeV/c for the 
K*. For higher pt the data tend to be flatter than the theory. 

The 0 meson with its mass 1.02 GeV/c^ is an important probe of fireball dynamics. 
If it is produced from the thermalized fireball at the common freeze-out, then its 
spectrum must be close to that of protons. In Fig. 12 we compare the data with 
predictions based on parameters from Table Data seem to be slightly steeper than 
the prediction. This is best seen in the ratios Ri plotted in the lower panels of the 
Figure. Such a disagreement might indicate slightly earlier freeze-out from the hreball. 
Note, however, that the deviation of the prediction from the data is much smaller than 
in the case of multistrange baryons and no indisputable sign of earlier decoupling can 
be claimed. 


5. Conclusions 

The blast-wave model with hadron production from resonance decays included can well 
reproduce single-particle spectra measured by ALIGE collaboration in Pb-|-Pb collisions 
at yJsNN = 2.76 TeV. Resonance contribution leads to shifts in the values of freeze-out 
parameters. For central collisions our freeze-out temperatures tend to be lower than 
those extracted by ALIGE m. On the other hand, the temperatures shift in the 
opposite direction for peripheral collisions. Shifts in {vt) have smaller relative size. 

Multistrange baryons show freeze-out at higher temperature and weaker transverse 
flow. This is consistent with earlier decoupling from the hreball. Note that such a 
behaviour is known from SPS [ 21 ] and RHIG [ 21 | data, as well. 

In hts to data from most peripheral collisions we obtain temperatures for kinetic 
freeze-out that are higher than what is usually assumed for the chemical freeze-out. On 
one side we can take this result as a parametrisation and the values of parameters are 
those which £t the data best. However, as soon as we want to interpret these values 
in the framework of some scenario for the fireball evolution, we may be witnessing 
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Figure 11. Comparison of measured K* spectra m with the predictions based on 
fitted parameters from Table Stars in the legend indicate that the centrality class 
of data 0-20% is compared to predictions for 0-5%, 5-10%, and 10-20%. To display 
all spectra in one figure we divide data from non-central collisions by factors 2, 4, 8, 
16, and 32. Lower panels show the ratios Rt = . 


limitations of the used model. 

In principle, our assumption of chemical equilibrium may seem inconsistent with 
individual normalisation of each spectrum as we performed it. If different spectra are 
multiplied with different factors this can be seen as effectively changing relative numbers 
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Figure 12. Comparison of measured ^ spectra m with the predictions based on 
fitted parameters from Table Stars in the legend indicate that the centrality classes 
20-30% and 30-40% of experimental data are compared to predictions based on the 
fit for 20-40% centrality in pions, kaons and (anti)protons. 
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of final state species. Caution would be particularly in place if we tried to predict 
the absolute yields of final hadron species. Here, however, we wanted to enhance the 
sensitivity to T and (vt) and thus we focussed on the shapes of the spectra which 
remain intact by the individual normalizations. Indeed, although rescaling one spectrum 
might have an influence on another one through changing yields of resonances which 
also contribute to a different hnal state species (like e.g. A contributing to pions and 
nucleons), this does not tie different spectra too strongly together. The reason is that 
the main contributing resonance decays are different for different hadron species: for 
pions the most important are p, u, rj [26]; for protons the A resonance; for A the 
and S* hyperons; for kaons the K*. Thus rescaling the yield for one species does not 
have much influence on another species. 

A more complete analysis would include also the model interpretation of the 
normalisations which are given by the volume. Such an analysis, however, must include 
femtoscopy and is planned for the future. 

Similar hts were performed recently with the Cracow single freeze-out model 
imiiH]. The main feature of that model is common chemical and thermal freeze- 
out with non-equilibrium fugacities introduced for strange and also non-strange species. 
Due to this, the freeze-out temperature htted there is generally higher than our results. 
Note however, that their freeze-out hypersurface has a different shape. Naturally, this 
opens the question if one can distinguish between the two models. Perhaps femtoscopy 
can provide the tool for this task. 

It is an intriguing speculation if the two approaches may be equivalent. Here we 
recall the scenario that matter hadronizes and cools down while keeping partial chemical 
equilibrium. This leads to the appearance of non-equilibrium chemical potentials for all 
species, i.e. chemical potentials that are not given by chemical potentials of conserved 
quantum numbers [2Z|. In the single freeze-out model, on the other hand, the system is 
born into the chemical non-equilibrium. We wonder how close to this non-equilibrium 
state the hreball comes in the cooling hreball scenario. 

A comprehensive picture of the gradual hreball cooling should also include earlier 
decoupling of multistrange baryons. 

Finally, let us come back to the issue that our results have been obtained in a 
scenario which assumes different chemical and thermal freeze-out and the difference 
between their temperatures is rather high. If one would revisit the question of the 
evolution between them, the inclusion of chemical potentials may be important. They 
may have inhuence on the shape of pt spectra especially for pions and we seem to observe 
the corresponding features in the data. This task will be addressed in a subsequent 
paper. 
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